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ABSTRACT 

We calculate the relative grain-grain motions arising from interstellar magnetohy- 
drodynamic (MHD) turbulence. The MHD turbulence includes both fluid motions and 
magnetic fluctuations. While the fluid motions accelerate grains through hydro-drag, 
the electromagnetic fluctuations accelerate grains through resonant interactions. We 
consider both incompressive (Alfven) and compressive (fast and slow) MHD modes and 
use descriptions of MHD turbulence obtained in Cho & Lazarian (2002). Calculations 
of grain relative motion are made for realistic grain charging and interstellar turbulence 
that is consistent with the velocity dispersions observed in diffuse gas, including cutoff 
of the turbulence from various damping processes. We show that fast modes domi- 
nate grain acceleration, and can drive grains to supersonic velocities. Grains are also 
scattered by gyroresonance interactions, but the scattering is less important than accel- 
eration for grains moving with sub-Alfvenic velocities. Since the grains are preferentially 
accelerated with large pitch angles, the supersonic grains will be aligned with long axes 
perpendicular to the magnetic field. We compare grain velocities arising from MHD 
turbulence with those arising from photoelectric emission, radiation pressure and H2 
thrust. We show that for typical interstellar conditions turbulence should prevent these 
mechanisms from segregating small and large grains. Finally, gyroresonant acceleration 
is bound to preaccelerate grains that are further accelerated in shocks. Grain-grain 
collisions in the shock may then contribute to the overabundance of refractory elements 
in the composition of galactic cosmic rays. 

Subject headings: ISM:dust, extinction — kinematics, dynamics — magnetic fields 



1. Introduction 

Dust is an important constituent of the interstellar medium (ISM) . It interferes with observa- 
tions in the optical range, but provides an insight to star-formation activity through far-infrared 
radiation. It also enables molecular hydrogen formation and traces the magnetic field via emission 
and extinction polarization (see reviews by Hildebrand et al. 2000, Lazarian 2000, 2002). The 
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basic properties of dust (extinction, polarization, etc.) strongly depend on its size distribution. 

The latter evolves as the result of grain collisions, whose frequency and consequences depend on 
grain relative velocities (see discussions in Draine 1985, Lazarian &; Yan 2002ab). 

Grain-grain collisions can have various outcomes, e.g., coagulation, cratering, shattering, and 
vaporization. For collisions with 5v < lO^^km/s, grains are likely to stick or coagulate, as the 
potential energy due to surface forces exceeds the initial center of mass kinetic energy. Coagulation 
is considered the mechanism to produce large grains in dark clouds and accretion disks. Collisions 
with Sv > 20km/s have sufficient energy to vaporize at least the smaller of the colliding grains 
(Draine 1985). It is likely that some features of the grain distribution, e.g., the cutoff at large 
size (e.g., Kim, Martin &; Hendry 1994), are the result of fragmentation (Biermann & Harwit 
1980). Even low- velocity grain collisions may have dramatic consequences by triggering grain 
mantle explosion (Greenberg &: Yencha 1973, Schutte & Greenberg 1991). 

Various processes can affect the velocities of dust grains. Radiation, ambipolar diffusion, 
and gravitational sedimentation all can bring about a dispersion in grain velocities. It is widely 
believed that, except in special circumstances (e.g., near a luminous young star, or in a shock wave), 
none of these processes can provide substantial random velocities so as to affect the interstellar 
grain population via collisions (Draine 1985), except for possibly enhancing the coagulation rate. 
Nevertheless, de Oliveira-Costa et al. (2002) speculated that starlight radiation can produce the 
segregation of different sized grains that was invoked to explain the imperfect correlation of the 
microwave and 100//m signals of the foreground emission (Mukherjee et al. 2001). If true it has 
important implications for the CMB foreground studies. However, the efficiency of this segregation 
depends on grain random velocities, which we study in this paper. 

The interstellar medium is magnetized and turbulent (see Arons Sz Max 1975, Scalo 1987, 
Lazarian 1999). Although turbulence has been invoked by a number of authors (see Kusaka et 
al. 1970, Volk et al. 1980, Draine 1985, Ossenkopf 1993, Weidenschilling & Ruzmaikina 1994) 
to provide substantial grain relative motions, the turbulence they discussed was not magnetized. 
Dust grains are charged, and their interactions with magnetized turbulence are very different from 
the hydrodynamic case. Lazarian &; Yan (2002a, henceforth LY02 and 2002b) applied the theory 
of Alfvenic turbulence (Goldreich Sz Sridhar 1995, henceforth GS95, see Cho, Lazarian & Vishniac 
2002a for a review) to grain acceleration and considered the motions that emerge due to incomplete 
coupling of grains and gas. Unlike the pure hydrodynamic case discussed by earlier authors, LY02 
took into account that the motions of grains are restricted by magnetic fields in the direction 
perpendicular to the field lines, and also took into account the anisotropy of Alfvenic turbulence. 

While Alfvenic turbulence is the turbulence in an incompressible fluid, the ISM is highly 
compressible. Compressible MHD turbulence has been studied recently (see review by Cho & 
Lazarian 2003a and references therein). In Yan &; Lazarian (2003, henceforth YL03) we identified 
a new mechanism of grain acceleration- gyroresonance -that is based on the direct interaction 
of charged grains with MHD turbulence. YL03 provided a test calculation of grain acceleration 
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in compressible MHD turbulence, both by hydro-drag and by gyroresonance (see also review by 
Lazarian & Yan 2003). 

In what follows, we describe grain acceleration by MHD turbulence in different phases of the 
ISM. Solving the Fokker-Planck equation including simultaneously the hydro-drag and gyroreso- 
nance would be a formidable task that we do not attempt here. Instead, wc try to simplify the 
problem by separating these two processes. For the random fluid drag, we use a simple scaling 
argument similar to the approach in Draine (1985) and LY02. While dealing with gyroresonance, 
we follow the approach adopted in YL03, i.e., we do not include the motion of ambient gas and 
assume that turbulence provides nothing but electromagnetic fluctuations. These approximations 
should yield correct answers when one of the mechanisms is dominant. When the accelerations 
arising from the two mechanisms are comparable the situation is more complicated, as the gaseous 
friction that we use for gyroresonance calculations will be, in general, affected by fluid motions. 
We do not develop explicit theory for this case. But taking into account that it is the motions at 
the decoupling scale that accelerate grains via hydro-drag, we think it is reasonable to estimate the 
velocity gains from the simultaneous action of hydro-drag and gyroresonance by adding them in 
quadrature. 

To describe the turbulence we use the statistics of Alfvenic modes obtained in Cho, Lazarian 
&; Vishniac (2002b, hereafter CLV02) and compressive modes obtained in Cho Sz Lazarian (2002, 
hereafter CL02, 2003bc)^ We apply our results to different phases of ISM, including the cold 
neutral medium (CNM), warm neutral medium (WNM), warm ionized medium (WIM), molecular 
cloud (MC) and dark cloud (DC) conditions, to estimate the implications for coagulation, shattering 
and segregation of grains. 

In what follows, we introduce the statistical description of MHD turbulence and damping 
processes (§2), describe motions arising from hydro-drag (§3) and gyroresonance (§4), apply our 
results to various ISM phases (§5), discuss astrophysical implications of our results (§6), and provide 
the summary in §7. 

2. MHD cascade and its damping 

MHD perturbations can be decomposed into Alfvenic, slow and fast modes (see Alfven &: 
Falthmmar 1963). Alfvenic turbulence is considered by many authors as the default model of 
interstellar magnetic turbulence. This is partially motivated by the fact that unlike compressive 
modes, the Alfvenic modes are essentially free of damping in a fully ionized medium (see Ginzburg 
1961, Kulsrud Sz Pearce 1969). Important questions arise. Can the MHD perturbations that charac- 
terize turbulence be separated into distinct modes? Can the linear modes be used for this purpose? 
The separation into Alfven and pseudo- Alfven modes is the cornerstone of the Goldreich-Sridhar 
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(1995, henceforth GS95) model of turbulence. This model and the legitimacy of the separation were 
tested successfully with numerical simulations (Cho & Vishniac 2000; Maron & Goldreich 2001; 
CLV02). Separation of MHD perturbations in compressible media into fast, slow and Alfven modes 
is discussed in GS95, Lithwick &; Goldreich 2001, and CL02. The actual decomposition of MHD 
turbulence into Alfven, slow and fast modes was performed in CL02, and Cho Sz Lazarian (2003, 
henceforth CL03) , who also quantified the intensity of the interaction between different modes (see 
below) . 

Unlike hydrodynamic turbulence, Alfvenic turbulence is anisotropic, with eddies elongated 
along the magnetic field (see Montgomery &; Turner 1981, Shebalin, Matthaeus, & Montgomery 
1983, Higdon 1984, Zank & Matthaeus 1992). This happens because it is easier to mix the magnetic 
field lines perpendicular to the direction of the magnetic field rather than to bend them. The 
GS95 model describes incompressible Alfvenic turbulence, which formally means that plasma (3 = 
-PgasZ-fmag; 1-6., the ratio of gas pressure to magnetic pressure, is infinity. The turbulent velocity 
spectrum is easily obtained. Calculations in CLV02 prove that motions perpendicular to magnetic 
field lines are essentially hydrodynamic. As the result, energy transfer rate due to those motions 
is a constant Ei^ ~ vf^jT}., where is the energy eddy turnover time ~ ("Ufe/^^)^^, where k±is 
the perpendicular component of the wave vector k. The mixing motions couple to the wave-like 
motions parallel to magnetic field giving a critical balance condition, i.e., k_\_Vk ~ ^yVA, where k\\ is 
the parallel component of the wave vector k, and Va is the Alfven speed. Prom these arguments, the 
scale dependent anisotropy fcy oc /c^and a Kolmogorov-like spectrum for the perpendicular motions 
Vk oc k~^/^ can be obtained. 

It was conjectured in Lithwick & Goldreich (2001) that the GS95 scaling should be approx- 
imately true for Alfven and slow modes in moderately compressible plasma. For magnetically 
dominated (i.e., low (3) plasma, CL02 showed that the coupling of Alfvenic and compressive modes 
is weak and that the Alfvenic and slow modes follow the GS95 spectrum. This is consistent with 
the analysis of HI velocity statistics (Lazarian &: Pogosyan 2000, Stanimirovic & Lazarian 2001) as 
well as with electron density statistics (see Armstrong, Rickett & Spangler 1995). Calculations in 
Cho & Lazarian (2003, hereafter CL03) demonstrated that fast modes are marginally affected by 
Alfven modes and follow acoustic cascade in both high and low (3 media. 

In what follows, we consider both Alfven modes and compressive modes and use the description 
of those modes obtained in CL02 and CL03 to study dust acceleration by MHD turbulence. 

The distribution of energy between compressive and incompressive modes depends, in general, 
on the driving of turbulence. CL02 and CL03 studied generation of compressive perturbations 
using random incompressive driving, obtaining an expression that relates the energy in fast ~ SV^ 
and Alfven ^ SV^ modes, 

{SVf/6VAf - 6VAVA/{Vi + C|), 
where Cs is the sound speed. This relation testifies that at large scales incompressive driving can 
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transfer an appreciable part of energy into fast modes. However, at smaller scales the drain of 
energy from Alfven to fast modes is marginal. Therefore the cascades evolve without much cross 
talk. A more systematic study of different types of driving is required. In what follows we assume 
that equal amounts of energy are transfered into fast and Alfven modes when driving is at large 
scales. 

We show that while simple scaling relations are sufficient for obtaining the velocities arising 
from hydro-drag, much more sophisticated tools are necessary for calculating gyroresonance (see 
Yan & Lazarian 2002, henceforth YL02, YL03). The corresponding statistics of turbulence is 
presented in Appendix B. 

At small scales the turbulence spectrum is truncated by damping. Various processes can 
damp the MHD motions. In partially ionized plasma, the ion-neutral collisions are the dominant 

damping process. In fully ionized plasma, there arc basically two kinds of damping: collisional or 
coUisionlcss damping (see Appendix A for details). Their relative importance depends on the mean 
free path / = vthT = 6 x 10^^{T/8000Kf ciR-^/n in the ISM (Braginskn 1965). If the wavelength is 
larger than the mean free path, viscous damping dominates. If, on the other hand, the wavelength is 
smaller than mean free path, then the plasma is in the collisionless regime and collisionless damping 
is dominant. 

To obtain the truncation scale, the damping time should be compared to the cascading time. 
As we mentioned earlier, the Alfvenic turbulence cascades over one eddy turn over time {k±Vk)~^ ~ 
(fey Va)~^. The cascade of the fast modes is a bit slower: 

Tk = u/k\l = (fc/L)-V2 X Vf/V', (1) 

where Vf is the phase velocity of fast waves, and V is the turbulence velocity at the injection scale 
(CL02) . If the damping is faster than the cascade, the turbulence is truncated. Otherwise, for the 
sake of simplicity, we ignore the damping and assume that the turbulence cascade is unaffected. 
According to CL02 the transfer of energy between Alfven, slow and fast modes of MHD turbulence 
is suppressed. This allows us to consider different components of MHD cascade independently. 

3. Grain Charge 

The net electrical charge on a grain is the result of competition between collisions with elec- 
trons, which add negative charge, and photoelectric emission and collisions with ions, which remove 
negative charge. We assume the grains to be spheres consisting of cither "astronomical silicate" or 
graphite, with absorption cross sections calculated as described by Wcingartner & Draine (2001a; 
henceforth WDOla), and photoelectric yields (as a function of Z) estimated by Weingartner & 
Draine (2001b; henceforth WDOlb). 

The grain charge depends on the electron density Ue- While many previous studies have 
assumed cosmic ray ionization rates ^ ~ 1 x 10~^^s~^ (e.g., Ruffle et al. 1998), recent observational 
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determinations (Black k, van Dishoeck 1991; Lepp 1992; McCall et al. 2003) suggest C « 1 — 
10 X 10~^^s~^ in H I clouds. We adopt an electron density ng ~ 0.03 cm~^ for CNM conditions, 
consistent with a detailed study of the ionization toward 23 Ori (Welty et al. 1999; Weingartner & 
Draine 2002), corresponding to an H ionization rate C, ~ 1.5 x 10^^^ s~^. 

For the outer regions of molecular clouds (MC) we take nn ~ 300 cm^'^ and ng/nn ~ 10"^, 
mainly due to photoionization of metals, where Gjjv = 0.1 is the UV intensity relative to the average 
interstellar radiation field. For "dark clouds" with nn ~ lO^cm"^, we consider Guv/ne = Icm^, 
resulting in negatively-charged grains. 

Fig. la shows the mean grain charge {Z) for graphite and silicate grains in various phases of 
the ISM (see Table 1). 

The charge on a given grain fluctuates. Let fz be the probability of the grain being in charge 
state Z, and let rz be the probability per unit time of leaving charge state Z. The characteristic 
time scale for the grain charge to fluctuate is tz = {{Z — (■^))^)/X^^ fzfz- Fig- lb compares the 

charge fluctuation time tz to the Larmor time (= Larmor period/27r), TLar = 1/^ = mg^c/ {Z)eB 
and the gas drag timescale t^^^^^ for subsonic motion. We see that, except for a < 10^® cm grains 
in dark clouds, the grain charge fluctuation time is much shorter than cither of these dynamical 
times, so that these fluctuations can be ignored and the charge on a given grain can be assumed to 
be constant, equal to the time-averaged charge {Z). 



4. Grain Motions arising from hydro-drag 

In hydrodynamic turbulence, the grain motions are caused by the frictional interaction with 
the gas. On large scales grains are coupled with the ambient gas, and the fluctuating gas motions 
mostly cause an overall advcction of the grains with the gas (Draine 1985). At small scales grains 
are decoupled. The largest velocity difference occurs on the largest scale where grains are still 
decoupled. Thus the characteristic velocity of a grain with respect to the gas corresponds to the 
velocity dispersion of the turbulence on the time scale ta- In the MHD case, the charged grains are 
subject to electromagnetic forces. If TLar > id, the grain does not feel the magnetic field. Otherwise, 
if TLar < ^di grain perpendicular motions are constrained by the magnetic field. 

As Alfvenic turbulence is anisotropic, it is convenient to consider separately grain motions 
parallel and perpendicular to the magnetic field. The motions perpendicular to the magnetic field 
are influenced by Alfven modes, while those parallel to the magnetic field are subjected to the mag- 
netosonic modes. According to §2, the perpendicular velocity fleld scales as ~ V {Tk/rmax)^^'^ > 
where Tmax = L/V is the time-scale on the injection scale. 

If the Larmor time TLar < Td^ grain perpendicular motions are constrained by the magnetic 
fleld. In this case, grains have a velocity dispersion determined by the turbulence eddy whose 
turnover period is ~ TLarj while grains move with the magnetic fleld on longer time scales. Since 
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Fig. 1. — left: \{Z)\ as a function of grain radius for carbon and silicate grains in 6 different 
environments: CNM,WNM,WIM,MC, DCl and DC2. \{Z)\ distribution is the same for DCl and 
DC2, which are referred to as DC in the plot, right: the gas-drag time t^ for subsonic grains, the 
Larmor time TLan and the charge relaxation time scale, tz, all as a function of grain size for silicate 
grains in the 6 different environments. 
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the turbulence velocity grows with the eddy size, the largest velocity difference occurs on the 
largest scale where grains are still decoupled. Thus, following the approach in Draine (1985), we 
can estimate the characteristic grain velocity relative to the fluid as the velocity of the eddy with 
a turnover time equal to TLar, 



= 0.9x 105cms-i(F5a5f/V(^^i0S^)'^'> (2) 

in which V5 = cms"-*^, 05 = a/10~^cm, Z = q/e, Liq = L/lOpc, = B/lfiG, and pgr is 

the mass density of grain. We adopt pgr = 3.8g cm~^ for silicate grains and pgr = 2gcm~^ for 
carbonaceous grains. 

Grain motions parallel to the magnetic field are induced by the compressive component of the 
slow mode with ~ 1^ x Th/Tmax^- For grain motions parallel to the magnetic field the Larmor 
precession is unimportant and the gas-grain coupling takes place on the translational drag time td- 
The drag time due to collisions with atoms t'^ = {apgr/nn){T^/8mnkBT)^^'^, where a is the grain 
size, nin is the mass of gas species, T is the temperature, is essentially the time for collision with 
the mass of gas equal to the mass of grain. The ion-grain cross-section due to long-range Coulomb 
force is larger than the atom-grain cross-section . For subsonic motions, the effective drag time 
^0 = t'Ja, where (Draine & Salpeter 1979) 



a 



1/2 



In 



3(fcBT)3/2 



aksT J 



2e3|Z|(7rxnH)^/2 



(3) 



where Xi is the abundance, relative to hydrogen, of ion i with mass rrii, x = Y^^Xi, and fz is the 
probability of the grain being in charge state Z. 

The characteristic velocity of grain motions along the magnetic field is approximately equal to 
the parallel turbulent velocity of eddies with turnover time equal to td 



v\\{a) = a ^—r \ . 



Pgr 



2tt 



"ninkBT 



= (3.8 X 10^ arv/s)a-^Via^/{nL^oTliQ), 
where Tioo = T/IOOK. Eq. (4) is valid for subsonic motion, v < {kBT/ninY^'^- 



(4) 



^We assume that turbulence is driven isotropically at the injection scale L. 
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When the grain velocity v relative to the gas becomes supersonic (Purcell 1969), the gas drag 
time td ~ ty{l + {9Tr/l28y/C^y/^. 

When TLar > td, grains are no longer tied to the magnetic field. Since at a given scale, the 
largest velocity dispersion is perpendicular to the magnetic field direction, the velocity gradient over 
the grain mean free path is maximal in the direction perpendicular to the magnetic field direction. 
The corresponding scaling is analogous to the hydrodynamic case, which was discussed in Draine 
(1985): 



v{a) = 



= a 



-1/2 
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( Pgr_ 
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(5) 



It is easy to see that the grain motions get modified when the damping time scale of the 
turbulence Tc is longer than either td or TXar- In this case, a grain samples only a part of the 
eddy before gaining the velocity of the ambient gas. In a turbulent medium, the shear rate dv/dl 
increases with the decrease of eddy size. Thus for Tc > max{t(i, TLar}, these smallest available eddies 
are the most important for grain acceleration. Consider first the perpendicular motions. If Vc is the 
velocity of the critically damped eddy, the distance traveled by the grain is A/ ~ X minjtd, 7T:,ar}- 
The shear rate dv/dl perpendicular to the magnetic field is r^^. Thus the grain experiences the 
velocity difi^erence in the direction perpendicular to the magnetic field 



v± Al X — — X mm{td,TLax}- (6) 
at Tc 

For the parallel motions, Al ^ Vc x td- From the critical balance in the GS95 model fcyVA ~ 
k±v± = T^^, the largest shear rate along the magnetic field should be dv/dl = fc^y ^ Vc/iVATc)- 
Therefore, in the parallel direction, the grain experiences a velocity difference Va/vc times smaller, 
i.e.. 



td 



The velocity dispersion induced by the compressional motion associated with the fast modes 
also causes motion relative to the ambient gas. The velocity fiuctuation for fast modes scales as 
V}^ oc A;^^/^ oc where u) is the frequency of fast modes. From similar considerations, we know 

that grains get velocity dispersions during minjrLar, idlji-c. , v ~ V(min{TLar5 ^dj/Tmaa;)^^^- Grains 
with min{rLar, ^cz} < Tc have reduced velocities v VcX ti^^^/tc ~ l^(rLar/Tmax)^''^(TLar/Tc)^/'', 
where Vc is the velocity of turbulence at the damping scale. Prom the scaling, we see that the 
decoupling from fast modes always brings larger velocity dispersions to grains than Alfven modes 
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{vk oc k~^^^) except for the situation when Alfven modes dominate MHD turbulence. The velocity 
fluctuations associated with fast modes are always in the direction perpendicular to B in low f3 
media (see Appendix D). Thus the grain velocities are also perpendicular to B. In high /? media, 
grains can have velocity dispersion in any direction as the velocity dispersions of fast modes are 
longitudinal, i.e., along k. 



5. Acceleration of Grains by Gyroresonance 

Gyroresonance acceleration of charged grains by a spectrum of MHD waves decomposed into 
incompressive Alfvenic, and compressive fast and slow modes (see CL02) was first described in 
YL02. The resonance happens when co — k^^vfi = nfi, (n = 0, ±1, ±2...), where lo is the wave 
frequency, fcy is the parallel component of wave vector k along the magnetic field, v is the grain 
velocity, fi is the cosine of the grain pitch angle relative to the magnetic field, and Q is the Larmor 
frequency of the grain. There are two main types of resonant interactions: gyroresonance acceler- 
ation and transit acceleration. Transit acceleration [n = 0) requires longitudinal motions and only 
operates with compressive modes. It happens when k^^vfj, = u, which requires particle speed to 
be super- Alfvenic v > Vf > Va- Although this condition is partially relieved owing to resonance 
broadening (see Yan &: Lazarian 2004) , transit acceleration of low speed grains is marginal because 
sub- Alfvenic particles can hardly catch up with the moving magnetic mirror. 

How can we understand grain gyroresonance? Gyroresonance occurs when the Doppler shifted 
frequency of the wave in the grain's guiding center rest frame LOgc = uj — k^^vfi is a multiple of 
the grain gyrofrequency. For low speed grains, we only need to consider the resonance at n = 1. 
The gyroresonance changes both the direction and absolute value of the grain's momentum (i.e., 
scatters and accelerates the grain). The efficiency of the two processes for charged grains can be 
described by the Fokker-Planck cocfRcicnts D^^ and Dpp/p^, where p is the grain momentum. The 
ratio of the two rates depends on the ratio of the grain velocity and the Alfven speed and pitch 
angle, p^D^^/Dpp = [{v{^/Va) + fj]'^, where C = 1 for Alfven waves and = fcy/A; for fast modes 
(see Appendix B). We see that the scattering is less efficient for sub- Alfvenic grains unless most 
grains move parallel to the magnetic field. We shall show later that as the result of acceleration, 
fi will tend to 0. Therefore in the zeroth order approximation, we ignore the effect of scattering 
and assume that the pitch angle cosine fi does not change while a grain is accelerated. In this case, 
the Fokker-Planck equation, which describes the diffusion of grains in momentum space, can be 
simplified (see Pryadko Petrosian 1997): 

dff^ dp ^ d . dp 

where / is the distribution function. Apart from acceleration, a grain is subjected to gaseous 
friction. Thus we describe the stochastic acceleration by the Brownian motion equation: 



- 11 - 



dv V 

where m is the grain mass, Y is the stochastic acceleration force, S = td/m is the mobihty coefficient. 
If we multiply equation(8) by v, and take the ensemble average, we obtain 

m^ = -Y + i^n. (9) 

The steady solution is achieved when the derivative on the left-hand is zero. Following an approach 
similar to that in Melrose (1980), we can get from Eq.(7) the energy gain rate for the grain with 
pitch angle 

r) = ^^(Vi?Pp(/^))- (10) 

The Fokker-Planck coefficient Dpp^jj) is calculated below. 

YL03 employed quasi- linear theory (QLT) to obtain Dpp{^) (see also YL02 and Appendix B) 
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{/ 



dk-^ 



Tu + {lO — kllVjJ, — O)^ 



,k±v± ,k±v± 



9, 
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(k) + e- 



£7^(k))]} 



(11) 



However, we should not integrate over all k because the contribution from large scales is spurious 
(see discussion in YL02). This contribution stems from the fact that in QLT, an unperturbed 
grain orbit is assumed, which results in non-conservation of the adiabatic invariant ^ = mv\/2Bo. 
Noticing that the adiabatic invariant is conserved when the electromagnetic field varies on a time 
scale longer than we truncate our integral range, namely, integrate from kres instead of 

the injection scale L~^. For Alfvenic turbulence to = |fc|||VA, the resonant scale corresponds to 
l^||,resl = ^/IVa — viJ,\. For fast modes, the resonant scale is kres = ^/|^/ — viJ,cos9\, where 
cos 9 = k\\lk. The upper limit of the integral k^ is set by the dissipation of the MHD turbulence, 
which varies with the medium. 

Integrating from fc^es to feci we obtain from Eq.(ll) and (10) the energy gain rate e as a function 
of V and /x. Then with e known, we estimate the grain acceleration. Solving Eq.(9) iteratively, we 



^Usually the real part is taken of the integral. However, we show in Appendix B that the integrand is real. 
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obtain the grain velocity as a function of time. We check that the grain velocities converge to 

a constant value after the drag time. As e increases with pitch angle, grains gain the maximum 
velocities perpendicular to the magnetic field and therefore the averaged fi decreases. This is 
understandable since the electric field, which accelerates the grain, is in the direction perpendicular 
to the magnetic field. 

6. Grain motions in the ISM 

Here we apply our results to various idealized phases of the interstellar medium. 

First consider a typical cold neutral medium (CNM), T = lOOK, nn = 30cm~^, rig = 
0.045cm~^, B = 6fiG, with corresponding va = km/s and f3 ^ 0.4. 

Our treatment of MHD turbulence requires that fluid velocities are smaller than the Alfven 
speed^. Therefore we assume that the injection of energy happens at an effective injection scale 
L where equipartition between magnetic and kinetic energies, i.e., V = Va, occurs. This effective 
injection scale can be different from the actual scale at which energy is injected. For instance, if 
we assume that the velocity dispersion at the scale I = lOpc is 5km/s, this means that turbulence 
in the CNM is super-Alfvcnic at this scale. The turbulence then follows a hydrodynamic cascade 
down to a scale L fs 0.64pc where the turbulent velocity becomes equal to Va = 2km/s; we identify 
this as the effective "injection scale" L^with injection velocity V = Va = 2km/s. Alternatively, it 
is possible that the turbulence at large scales proceeds in tenuous warm media with Alfven speed 
larger or equal to 5km/s. Nevertheless, the statistics of fast modes will not be changed in the CNM. 

In partially ionized media, the damping is dominated by the viscosity arising from neutrals. 
The turbulence is assumed damped^ when its cascading time scale = tdamp'i this defines the 
cutoff scale /cy ^ = 4 x 10~^^cm~^ for Alfven modes and kc = 7 x 10~^^cm'~-'^ for fast modes (see 
Appendix A). Assuming that the grain velocities are smaller than the phase speed of fast modes, 
we find that the prerequisite for gyroresonance kc > kres is the same as TLar > Tc-, the condition for 
effective hydro drag (see Fig.l). For a silicate grain, the critical grain size Oc ~ 4 x 10~^cm for fast 
modes and ac ~ lO^^cm for Alfven modes. Grains smaller than the critical size are not effectively 
accelerated by the corresponding turbulent mode. 

The CNM is a low /3 medium, so the correlation tensors for the low /3 case are applied. As has 



''Otherwise magnetic field is not dynamically important. Turbulence is essentially hydrodynamic. See the following 
discussions. 

^This picture is not self-consistent as we expect to have turbulent generation of magnetic field which will bring 
kinetic and magnetic energy to equipartition at the injection scale (see arguments in Cho, Lazarian & Vishniac 
2002a), but the idea of super- Alfvenic turbulence percolates in the literature (Boldyrev, Nordlund & Padoan 2002). 

^Thus wG ignore the effect of slowly evolving magnetic structures associated with a recent reported new regime of 
turbulence below the viscous damping cutoff (Cho, Lazarian, & Vishniac 2002c, Lazarian, Vishniac, & Cho 2004). 
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been discussed in YL02, the interactions with Alfven modes are less efficient than those with fast 
modes because of the anisotropy of Alfven modes. Thus we shall consider only fast modes for later 
calculations. The gyrorcsonance with fast modes in the CNM can accelerate grains to supersonic 
velocities, e'^ increases with pitch angle. If we average e over fi, we will get mean velocities which 
are smaller than the maximum values by less than 20%. In Fig.l we plot the velocity of grains with 
pitch angle equal to 90° as a function of grain size since all the mechanisms preferentially accelerate 
grains in this direction. 

For the warm neutral medium (WNM), T = 6000K, nn = O.Scm^^, Ue = O-OScm^^, B = 
5.8/xG. Assume the velocity dispersion V = Va = 20km/s at the injection scale L = lOOpc. 
Turbulence is mainly subjected to neutral-ion damping. The fast modes are cut off at /cc = 4 x 
10"^'' cm~^ (see Appendix A). Comparing kc with kres, we find Oc ^ 2 x lO^^cm for silicate grains. 
The WNM has (3 ~ 0.25, so we use the tensor given in Eq.(B4) for fast modes. Integrating from 
kres to kc and solving Eq.(9), we obtain grain velocities. The maximum values are shown in Fig.2. 
We see large grains can be accelerated to supersonic speeds^. The fact that these grains approach 
the Alfven speed makes our approximation that acceleration dominates scattering less accurate, 
but the result is correct within a factor of unity. Smaller grains are accelerated only by the hydro 
drag, which is far less effective. 

The warm ionized medium (WIM) has T = 8000K, rig = 0.1 cm~^, B = 3.35^G, with corre- 
sponding f3 ~ 0.33. The injection scale and speed are the same as in the WNM. The WIM is fully 
ionized and in low /3 regime. Fast modes, in this case, are mainly affected by collisional damping. 
This damping increases with 6,^ and doesn't exist for parallel modes (see Appendix A). Thus there 
are always modes interacting with grains though the energy available is less at smaller scales. Fol- 
lowing the same routine as above, we get the grain velocities. We see from Fig.4 the nonmonotonic 
variation of grain velocity with the size. This arises from the fact that the charging for grains in 
the WIM has a complex dependence on grain size (see also Fig. la). 

Molecular cloud (MC) gas has T = 25K, nn = 300cm~'', and we adopt a magnetic field 
strength B = ll^uG as suggested by observations (Crutcher 1999), corresponding to Alfven speed 
Va = 1.2 kms ~^ . The injection scale is taken to be L = Ipc and the injection velocity is V = Va- 
The damping scale (see Appendix A) of the turbulence is kc = 4.5 x 10~^^ cm corresponding to 
resonant scales of silicate grains with a = 8 x 10~^cm. By following the same procedure, we obtain 
the grain velocity distribution as shown in Fig.5. 

We consider a typical dark cloud (DC) with T = lOK, nn = lO^cm"^, and B ~ 80/iG, 
corresponding to Va = 1.5kms~^. The injection scale L = Ipc and velocity V = Va- The 



^Unlike hydro-drag the gyroresonanco can potontiaUy accelerate grains to velocities much higher than the velocity 
of turbulent motions. For typical ISM conditions this, however, does not happen. 

®This 6 dependence makes the treatment of damping more complicated if taking into account field line wandering 
(see Yan & Lazarian 2004). 
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ne(cm~^) 


0.03 


0.03 
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0.03 
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0.001 


Guv 


1 


1 


1 


0.1 


0.01 


0.001 


B(^G) 


6 


5.8 


3.35 


11 


80 


L(pc) 


0.64* 


100 


100 


1 


1 


y = Vk(km/s) 


2* 


20 


20 


1.2 


1.5 


drunpiiig 


iieulrMl-it)u 


iieuLrMl-iou 


c;t)llisi(niMl 


iieulrMl-iou 


icni-iKuU ral dccDupliug 


A;c(cin"^) 


7 X 10-^5 


4 X 10"^^ 


NA 


4.5 X 10"^^ 


5.3 X 10"^^ 


5.3 X 10"^^ 



Table 1: The parameters of idealized ISM phases and relevant damping. Among them, nn is 
the number density of H, Ue is the number density of electron, Guv is the UV intensity scale 
factor, L is the injection scale of fast modes, V is the injection velocity. The dominant damping 
mechanisms for fast modes are given with the corresponding damping timescale Tc- CNM=cold neu- 
tral medium, WNM=warm neutral medium, WIM=warm ionized medium, MC=molecular cloud, 
DC=dark cloud. * See text for the explanation of smaller L and V for CNM. 
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Fig. 2. — Relative velocities as a function of radii (solid line) in CNM, left: for silicate grains, right: 

for graphite grains. The dotted lines represent the gyroresonance with fast modes. Gyroresonance 
works only for large grains owing to the cutoff by viscous damping. The cutoff scales for fast and 
Alfven modes are different due to their different scalings and the anisotropy of Alfven modes. The 
dashed lines are the result from hydro drag by Alfven modes (see LY02) , the dashdot lines represent 
the hydro drag by fast modes. Contributions from different mechanisms are approximately additive 
in squares, i.e., v^^^ = T^ivf (heavy solid lines). The grain velocity driven by H2 formation (dashdot 
line) is plotted to illustrate the issue of grain segregation in CNM (see text). The part marked by 
open circles is nonphysical because thermal flipping is not taken into account. 
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Fig. 4. — Same as Fig.2, but in the WIM. The oscillations in these curves are due to the variation 
in charging of grains. 
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ionization in DC is so low that the fluid becomes decoupled in the middle of the cascade, where 
the ion-neutral collision rate per neutral is equal to the turbulence decay rate t^^- Below this 
decoupling scale, neutrals will not follow ions and magnetic fields, and the turbulence becomes 
hydrodynamic. In view of the uncertainty of the cosmic-ray ionization rate, we adopt two models 
DCl and DC2 with electron densities rig = 0.01 cm"^ and ng = 0.001 cm~^. The grain charge 
distribution is the same for both DCl and DC2, because we assume the same Guv/ne- The main 
difference is the decoupling scale of the MHD cascade. Combining Eq.(l) and Eq.(A2), wc can 
obtain the decoupling scales kc = 5.3 x 10~^^cm~^ for DCl and kc = 5.3 x 10~^''cm~^ for DC2, 
corresponding to silicate grain size ~ 3 x lO^^cm for DCl and Cc ~ 2 x 10~^ cm for DC2. By 
following the same procedure, we obtain the grain velocity distribution as shown in Fig. 6 and 
Fig. 7. 

It is shown in Fig. 5,6&:7 that the acceleration by gyroresonance in both MC and DC is not 
as effective as in the lower density media, for two reasons. First, the low levels of UV and low 
temperatures result in reduced grain charge (see Fig. la). Secondly, because of the increased 
density, the drag time td is reduced. 

The grain velocities in CNM found here are smaller than in an earlier calculation (see YL03) 
because we adopt smaller values for magnetic field B and injection velocity V. 

It should be noted that the strength of magnetic fields in the ISM is still somewhat uncertain 
and may vary from place to place. We adopted a particular set of values in above calculations. 

How would the results vary as the magnetic field strength varies? First of all, we know that the 
critical condition for acceleration is kres > k^: grains with kres < k^ cannot be accelerated. Thus 
the cutoff grain radius Uc varies with the medium, Oc ~ [3q{ac)TcB /Att pj^^^ ; grains with a < Uc are 
not subject to gyroresonant acceleration. 

The magnitude of the velocity is a complex function of the magnetic field. For illustration, in 
Fig. 7 we show the grain velocities calculated for magnetic fields a factor of 3 stronger or weaker 
than the values in Table 1 for the CNM and DCl environment. Since the hydro drag by fast 
modes decreases with the magnetic field, the relative importance of gyroresonance and hydro drag 
depends on the magnitude of the magnetic field. In magnetically dominant regions, gyroresonance 
is dominant. In weakly magnetized regions, the frictional drag provides the highest acceleration 
rate. The injection scale is another uncertain parameter, but the grain velocity is not so sensitive 
to it provided that the injection scale is much larger than the damping scale. 



Discussion 



Shattering and Coagulation 



With the grain relative velocities known, wc can make predictions for grain shattering and 
coagulation. For shattering, we adopt the Jones et al. (1996) results, namely, for equal-sized 
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Fig. 6. — Same as Fig.2, but in DCl. 
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Fig. 8. — Relative velocities gained from gyroresonance as a function of radii for different magnetic 
field strength, left: in CNM, right: in DCl. Solid lines are the results for the precedent values of 
magnetic field. Dashed lines refer to the results with 3 times stronger magnetic field. Dash-dot 
lines represent the cases with 3 times weaker magnetic field. 
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particles the shattering threshold is 2.7km/s for siHcate grains and 1.2km/s for carbonaceous grains. 
The critical sticking velocity was given by Chokshi et al. (1993) (see also Dominik & Tielens 1997)^: 



2.14F,tick 



4 + 4 1'/' 7^/^ 



(ai+a2)3j EV3i25/6^1/2' 



is the maximum relative velocity for coagulation of equal-size spherical grains, where 7 is the surface 
energy per unit area, R = 0102/(01 +02) is the reduced radius of the grains and E is related to 
Poisson's ratio Ui and Young's modulus Ei by l/E = [(1 — vif /Ei + (1 — ^2)^ / E2] and we have 
introduced a factor FgUck ~ 10 since the experimental work by Blum (2000) shows that the critical 
velocity is an order of magnitude higher than the theoretical estimate of Chokshi et al. We use 7, 
E and u for Si02 and graphite from Table 1 of Dominik &; Tielens (1997), and consider collisions 
between equal-size grains (oi ~ 02). Comparing these critical velocities with the velocity curve we 
obtained for various media, we can get the corresponding critical size for each of them (see Table 
2). 

Correlation between turbulence and grain sizes 

The grain velocities are strongly dependent on the maximal velocity of turbulence V at the 
injection scale, which is highly uncertain. The critical coagulation and shattering sizes thus also 
depend on the amplitude of the turbulence. Variations in the level of turbulence could lead to 
regional differences in the grain size distribution. 

Elements in cosmic rays 

It has been shown that the composition of galactic cosmic rays appears to be correlated with 
elemental volatility (Ellison, Drury &; Meyer 1997). The more refractory elements are systematically 
overabundant relative to the more volatile ones. This suggests that the material locked in grains 
must be accelerated more efficiently than gas-phase ions (Epstein 1980; Ellison, Drury &: Meyer 
1997). The stochastic acceleration of grains, in this case, can act as a preacceleration mechanism. 
The ions released from the grains in the shock by sputtering or in grain-grain collisions can then be 
further accelerated in the shock, and explain the overabundance of refractory elements in galactic 
cosmic rays. 

Heavy Element Depletion and Grain Alignment 



®Note a misprint in the exponent of Young's modulus in eq.(lO) of Dominik & Tielens (1997). 
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Table 2: Size ranges for shattering and coagulation in different medium. NA=not applicable. 
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Our results indicate that grains can become supersonic through interaction with fast modes. 
Grains moving with velocities larger that those of heavy ions could sweep up heavy elements which 
may be advantageous from the point of explaining observations (Wakker & Mathis 2000). Our 
calculations show that while such velocities are readily achievable, the sign of charging may present 
a problem for such "vacuum cleaning" of the ISM. For instance, silicate grains in MC can be 
accelerated to > 4 x 10^ cm/s, which is larger than the thermal speed of heavy ion. Therefore 
the capture rate for ions by positively charged grains (< 2 x 10^^ cm) would be increased. Grains 
smaller than 2 x 10~^ cm will be negatively charged in the MC. For such grains the cross section 
for Coulomb capture of ions will decrease. If such small grains retain captured ions, this would 
result in a decrease of the rate of depletion of metals on grains. The actual rates of depletion on 
fast moving grains are important and will be identified elsewhere for particular phases of the ISM. 

Grains moving supersonically can be aligned mechanically (see a review by Lazarian 2003 and 
references therein). As pointed out earlier, the scattering is not efficient for slowly moving grains 
so that we may ignore the effect of scattering on the pitch angle. Since the acceleration of grains 
increases with the pitch angle of the grain(scc Eq.(lO) and (11)), the supersonic grains will tend 
to have large pitch angles. As first discussed by Gold (1952), gas drag acting on these grains will 
tend to cause them to spin with angular momenta perpendicular to their motion, and therefore 
tending to be parallel to the magnetic field direction. Dissipational processes will tend to orient the 
spinning grains with their long axes perpendicular to their angular momentum, resulting in grain 
alignment with long axes perpendicular to the magnetic direction. 

Grain Segregation and Turbulent Mixing 

Our results are also relevant to grain segregation. Grains are the major carrier of heavy 
elements in the ISM. The issue of grain segregation may have significant influence on the ISM 
metalicity. Subjected to external forcing (WDOl, Ciolek and Mouschovias 1996), grains gain size- 
dependent velocities with respect to gas. WDOl considered the forces on dust grains exposed 
to anisotropic interstellar radiation fields, including photoelectric emission, photodesorption, and 
radiation pressure, and calculated the drift velocity for grains of difi'erent sizes. The velocities 
they got for silicate grains in the CNM range from O.lcm/s to lO^cm/s. Grains can move along 
magnetic field lines due to the uncompensated forces, e.g. due to active sites of H2 formation 
(see P79; Lazarian &; Yan 2002b)^° Fig. 2a shows that the turbulence produces larger velocity 
dispersions^^ . Those velocities are preferentially perpendicular to magnetic field, but in many cases 
the dispersion of velocities parallel to magnetic field will be comparable to the regular velocities 
above. This dispersion stems from both the fact that the transpositions of matter by fast modes 
arc not exactly perpendicular to magnetic field (see plot in Lazarian & Yan 2002b) and due to 



^"These forces would be mitigated in molecular clouds, which would induce inflow of dust into molecular cloud. 
The latter would affect metalicity of the newborn stars. 

^^Our calculation show that for the chosen set of parameters the effects of systematic thrust are also limited (see 
LY02, Lazarian & Yan 2002b). 
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randomization of directions of grain velocities by magnetized turbulence (YL03). 

More important is that if reconnection in turbulent medium is fast (see Lazarian & Vishniac 
1999, Lazarian et al. 2004), the mixing of grains over large scales is provided by turbulent diffusivity 
~ VL/4:. Usually it was assumed that the magnetic fields strongly suppress the diffusion of charged 
species perpendicular to their directions. However, this assumption is questionable if we notice that 
motions perpendicular to the local magnetic field are hydrodynamic to high order as suggested by 
Cho, Lazarian &; Vishniac (2002). Recent work by Cho et al. (2003) found that the diffusion 
processes in MHD turbulence are almost as effective as in the hydrodynamic case if the mean 
magnetic field is weak or moderately strong (i.e., B < the equipartition value) which would imply 
that grains can be mixed by the MHD turbulence. Lazarian Sz Yan (2003) therefore concluded 
that the segregation of very small and large grains speculated in de Oliveira-Costa et al. (2002) is 
unlikely to happen for typical interstellar conditions. 

8. Summary 

We have calculated the relative motions of dust grains in a magnetized turbulent fluid. It has 
been known for decades that turbulence can give rise to significant grain-grain velocities. However, 
earlier treatments disregarded the magnetic field and used Kolmogorov turbulence. Magnetohydro- 
dynamic (MHD) turbulence includes both fluid motions and magnetic fluctuations. While the fluid 
motions bring about decoupled motions to grains, the electromagnetic fluctuations can accelerate 
grains through resonant interactions. 

Calculations of grain relative motion are made for different phases of the ISM with realistic 
grain charging, and with turbulence that is consistent with the velocity dispersions observed in 
interstellar gas. We account for the cutoff of the turbulence from various damping processes. We 
show that fast modes dominate grain acceleration, and can drive grains to supersonic velocities. 
Grains are also scattered by gyroresonance interactions. The scattering rate is less efficient than 
acceleration for grains moving with sub-Alfvenic velocities. 

Since the grains are preferentially accelerated with large pitch angles, the supersonic grains 
tend to be aligned with long axes perpendicular to the magnetic field. 

Gyroresonant acceleration is bound to preaccelerate grains that will then be further accelerated 
by shocks. Grain-grain collisions and sputtering in the shocks will inject suprathermal ions which 

can then undergo further acceleration in the shock, potentially accounting for the observed excess 
of refractory elements in the composition of galactic cosmic rays (e.g., Epstein 1980; Ellison et al. 
1997) 
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from NSF grant AST0243156 and that from the Center for Magnetic Self-Organization in the 



- 22 - 



Astrophysical and Laboratory Plasmas. BTD acknowledges partial support from NSF grant AST- 
9988126. 



A. Damping of MHD waves 

Below we summarize the damping processes that we consider in the paper. 
Neutral-ion damping 

In partially ionized medium, a combination of neutral viscosity and ion-neutral collisional 
coupling provides damping (see LY02). If the mean free path for a neutral, Z„, in a partially ionized 
gas with density ntot = nn + rii is much less than the size of the eddies under consideration, i.e. 
Ink <^ 1, the damping time 

tdamp ~ K'k-^ ~ ( — ) ilnVn)-'k-\ (Al) 

where is effective viscosity produced by neutrals^^, is the thermal velocity of the neutrals, 
and the mean free path of a neutral Z„ is influenced both by collisions with neutrals and with 
ions. The rate at which neutrals collide with ions is proportional to the density of ions, while the 
rate at which neutrals collide with other neutrals is proportional to the density of neutrals. The 
momentum transfer rate coefficient for neutral-neutral collisions is ~ 1.7 x 10^"'^'^(T/K)*'"^ cm'^ s^^ 
(Spitzer 1978), while for neutral-ion collisions it is ~ (fr^m) ~ 1.9 x 10~^ cm^ s~^ (Draine, Roberge 
h Dalgarno 1983). Thus collisions with other neutrals dominate for ni/rin less than ~ 0.09r°-^. 

Effects of charged grains 

Magnetic perturbations can get decoupled from the fluid motions because neutrals are imper- 
fectly coupled to the ions in partially ionized medium^^. The coupling between ions and neutrals 
is determined by the ion-neutral collisional rate: 

*n/ = -7 M^rf^in) (A2) 

where Vr is the ion-neutral relative velocity, ain is the ion-neutral collisional cross section, rrii and 
m„ are the typical ion and neutral masses, is the ion number density. When the collisional 
time tni is equal to the wave period, neutrals are decoupled from magnetic field, and turbulence 
becomes hydrodynamic. In molecular clouds, grains can take substantial portion of the total charge. 
The contribution of charged grains to coupling neutrals to magnetic fields depends on grain size 



The viscosity due to ion-ion collisions is typically small as ion motions are constrained by the magnetic field. 

^^We do not discuss here a viscosity-damped regime of MHD turbulence that takes place in the partially ionized 
gas below the scale at which viscosity damps kinetic motions associated with magnetic field (see theory of this regime 
in Lazarian, Vishniac & Cho 2004) 
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spectrum. The ratio of the ion-neutral coUisional rate to the grain-neutral collisional rate is (Nishi 
& Nakano 1991, Elmegreen & Fiebig 1993): 



tj ^[o.25n//^ (MRN) .^3. 
{ng(Tg)vn \3.3n4^/^ (MW) ' 

where MRN refers to the grain size distribution proposed by Mathis et al. (1977) and MW stands for 
the distribution suggested by Mathis & Whiffen (1989). Prom this expression, we can see that ions 
are always the dominant contribution for the coupling in MC. In DC, the situation will depend on 
grain size distribution. DC is a denser region where observations favor MW distribution (Elmegreen 
& Ficbig 1993). Thus presumably the contribution from grains in DC is also subdominant and we 
neglect it in the main text. 

Collisionless damping 

The nature of collisionless damping is closely related to the radiation of charged particles in 
magnetic field. Since the charged particles can emit plasma waves through acceleration (cyclotron 
radiation) and Cherenkov effect, they also absorb the radiation under the same condition (Ginzburg 
1961). The damping rate 7d = 'T'd^ of the fast modes of frequency u; for /3 <C 1 and ^ ~ 1 (Ginzburg 
1961) is 



4 cosp 



^^P R 2^ +^®^P(~^5 YB' 

mu ITiYlp cos^ & p cos^ a 



(A4) 



where nie is the electron mass. The exact expression for the damping of fast waves at small 9 was 
obtained in Stepanov^^ (1958) 
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where is the ion gyrofrequency. 
Ion viscosity 



2^li/(3 for k > ^i/fiVA 



(A5) 



In a strong magnetic field (OjTj ^ 1) the transport of transverse momentum is prohibited 
by the magnetic field (along z). Thus transverse viscosity r/j, is much smaller than longitudinal 



*We corrected a typo in the corresponding expression. 



viscosity 770, rf± ~ "Ho / {^iTi)^ ■ Following Braginskii (1965), we can find the damping rate is(see Yan 
& Lazarian 2004): 



r _\k'im/^Pi for/3<l 

U2r,o(l-3cos2e)76p, for^»l ^ ' 



For more discussion, see Yan &: Lazarian (2004). 



B. Fokker-Planck coefficients 

In quasi-linear theory (QLT), the effect of MHD waves is studied by calculating the first order 
corrections to the particle orbit in the uniform magnetic field, and the ensemble-averaging over the 
statistical properties of the MHD waves (Jokipii 1966, Schlickeiser & Miller 1998). Obtained by 
applying the QLT to the collisionless Boltzmann-Vlasov equation, the Fokker-Planck equation is 
generally used to describe the evolvement of the gyrophase-averaged particle distribution. 
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where p is the particle momentum. The Fokker-Planck coefficients D^^,D^p,Dpp are the funda- 
mental physical parameter for measuring the stochastic interactions. 
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fcmax corresponds to the dissipation scale, 7?., L refer to left- 



and right-circularly polarized modes, and ^ = tan ^ kx/ky. 
The correlation tensors are defined as following: 

< S„(k, t)S^(k', t + T)> /Bl = <)(k - k')M«;3(k)e-^/^'= 

< v^{^^,t)B*p{^^,t + t) > /VaBo = <5(k - k')Cc,^(k)e-^/^'= 

< Va{Kt)v*^{^',t + T) > IVl = 6ik - k')ifa/3(k)e-^/^S 



(B2) 
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where -Bq,^^, w^^^ are respectively the magnetic and velocity perturbation associated with the 

turbulence, Tfc is the nonlinear dccorrclation time and essentially the cascading time of the turbu- 
lence. For the balanced cascade wc consider (sec discussion of our imbalanced cascade in CLV02), 
i.e., equal intensity of forward and backward waves, Cij(k) = 0. 

The magnetic correlation tensor for Alfvenic turbulence is (CLV02), 
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(B3) 



where lij = {5ij — kikj/k"^} is a 2D tensor in ,x — y plane which is perpendicular to the magnetic 
field, L is the injection scale, V is the velocity at the injection scale. Slow modes are passive and 
similar to Alfven modes. The normalization constant is obtained by assuming equipartition = 
/ dk^ Si=i MuBq/Stt ~ Bq/Stt. The normalization for the following tensors below are obtained in 
the same way. 

According to CL02, fast modes are isotropic and have one dimensional energy spectrum E{k) oc 
k~^/^. In low P medium, the corresponding correlation is (YL03) 
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where 9 is the angle between k and B, Hij = kikj/k\ is also a 2D tensor in a; — y plane. The 
factor cos^ 6 represents the projection as magnetic perturbation is perpendicular to k. This tensor 
is different from that in Schlickeiser & Miller (1998). For isotropic turbulence, the tensor of the 
form oc Ek{Sij — kikj/k"^) was obtained to satisfy the divergence free condition k • JB = (see 
Schlickeiser 2002). Nevertheless, the fact that (5B in fast modes is in the k-B plane places another 
constraint on the tensor so that the term 5ij doesn't exist. 

In high (3 medium, fast modes in this regime are essentially sound waves compressing the 
magnetic field (Goldreich &; Sridhar 1995, Lithwick & Goldreich 2001, CL03). The compression of 
magnetic field depends on plasma The corresponding x-y components of the tensors are 



Tfe = (A; X kmvn)-^'^ X Cs/V\ (B5) 

where Cg is the sound speed. The velocity perturbations in high (3 medium are longitudinal, i.e., 
along k, thus we have the factor sin^ 6 and also a factor V^/C^ = 2//3 from the magnetic frozen 



A%(k) 
i^,,(k) 



,-7/2 COS^e//? 

^^sm eH,,k 



condition ujS'B ~ k x (v^ x B). We use these statistics to calculate grain acceleration arising from 
MHD turbulence. 

The spherical components of the correlation tensors are obtained below. 
For Alfven modes, their tensors are proportional to 



/ sin^^ -cos^sin^ \ 
\ — COS (j) sm (f) cos / 



Thus we get 



and 



For fast modes, their tensors have such a component 



inn - J^cc - — — - -^{.^xx + J-yy) - 2' 



'^Icn = X e^2^ + (^L±!^ X e^^^ = (4. - cos 20 + (4, + sin 20 



H,,=Ak-^^ ^°f'^^ cos0sin0y 
\ cos sm sm / 

Thus we have 



H-nn = Hcc = -^^H^x + Hyy) = -, 

and 



(^''^Hnc + e-'^'f'Hcn = (H^x - Hyy) cos 20 + {H^y + Hy^) sin 20 = 1. 
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C. Drag due to the dipole moment of grain 

The plasma drag to the rotational motion for stationary grains has been considered earlier 
(Anderson &; Watson 1993, Draine & Lazarian 1998, henceforth DL98). Similarly, for a grain with 
electric dipole moment /x, there also exist forces between the grain and nearby ions. Consider the 
effects in the comoving frame of grain. In this frame ions move at speed v with impact parameter 
b. To simplify the problem, we define a stopping cross section as = irbQ, where 6o is defined by 
Zie/i/bQ = rriiV^. For impact parameters b < max(a, 6o)) the interaction is strong. We will assume 

a<b,= [^\'' = 5.4 X 10-^ cm Z^/^i ( ^) ^ (CI) 



where jii = n/lO debye. If we assume the ions to be scattered isotropically, then the drag force 
due to strong scattering events is 

F,^n,Z,e„^ ;f (C2) 

6 max[uTO,(4/3V7r)VprJ 

where Vm = (2fcT/m)^/^ is the most probable speed of ions at temperature T. For ions with impact 
parameter > bo, we assume their trajectories are barely changed during the collisions with the 
grain. Define the direction of v as the polar axis e^, let pericenter be at (6,0,0), and let t be the 
time from pericenter. The force on the ion from the dipole moment is 

= J^^^^^i^^l^^C^^^ + ^^^^ - ^^^^) + ^zl^^i'^vH'^ + 3vtb - b^) - eyiiy{b^ + vH'^)). (C3) 

Integrated over time from — oo to oo, this expression yields the total momentum delivered to the 
grain: 

Ap = -^(-e^Mx + ej/M2/)- (C4) 

p increases in a random walk fashion, therefore the impulses of individual collisions should be 
added in quadrature. If we now average over random orientation of ji and then integrate over 
impact parameters and thermal distribution of ion velocities, we find 



dvAnv'^ fi{v)v / 2'Kbdb- ( ,„ ^ ) = — - — niZieiimiVjn , (C5) 
Jbo 3 \ b-'v J 3 



where for the upper cutoff we take the Debye length 62 = (kT/ATrriee'^)^/'^ ^ 6o- Using the 
fluctuation-dissipation theorem, we then can get the damping force for subsonic grains, 

F = F, + ^^ = ^-^n..Z,e/-^ . (C6) 
6kT dt 9 ^Vm 



^Hn DL98 it was estimated that /ii « 0.93(o/10"'' cm)*/^. 
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If the grains becomes supersonic, the fluctuation-dissipation theorem is no longer apphcable. 
Since the interaction is elastic, the loss of the momentum in the direction parallel to the moving 
direction can be obtained by assuming energy conservation: p^+Pi. ~ const. The dipole interaction 
is weak interaction so that Ap/p <C 1 during one encounter. Thus we can estimate the momentum 
loss in the z direction as Ap^ = Ap^/2p, where Ap±_ is given by Eq.(C4). Then integrating over 
impact parameter, we get the damping force 

nC^'^ 4 ( Z eu\ ^1 Ttt 

F = TTHiZiefi + / 27rM6 - — = —rnZiefi . (C7) 

To determine the importance of this force, we compare it with other drag forces. Using the 
dipole moment estimated by DL98, we find the dipole drag is smaller than collisional drag. However, 
for very small neutral grains with a dipole moment in ionized gas, dipole plasma drag may play a 
more important role. 

We note that the estimate by DL98 of rotational excitation by "plasma drag" acting on the 
grain dipole moment overestimated the transfer of angular momentum by using the weak interaction 
approximation for all trajectories with b > a. Assuming random scattering, the mean square angular 
momentum transfer from strong scattering events will be ~ 2{mivb)^, and thus the contribution of 
impact parameters a < 6 < 6o to dL? /dt is 



dl? 
dt 



oo 

rii I dvATTv'^fj 
'o 

Zfe'„' 



/•bo 

{v)v / : 

J a 



1 



2Trbdb 2{mvby = Am^Uim^ \ dvv^fi{v){b^ - a*) 



where T2 = T/IOOK, a_7 = a/10 ^cm. For comparison, for impact parameters a < 5 < 60 DL98 
found 



dl? 
dt 



32V^ zy^i^ 



3 v„ 



bo 



In ( — 

a 



1/2 



4.07 



1/4^1/2 



(C9) 



This part of contribution is comparable with the total if a_7 < 1 and T2 < 1 or if the angle 
between the dipole moment and the rotation velocity is close to 90° (see eq.(B35) in DL98). For 
instance, at T = lOOK, for grains with a_7 = 0.5, the part given by eq.(C9) is ~ 35% of the total 
for cos^ = 1/3. In such cases, the correction owing to the strong scattering as given in Eq.(C8) 
should be taken into account. 
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D. Angle between B and v 

A fundamental question arises from the fact that in MHD turbulence wave vectors are not 
aligned along magnetic field lines, as is the case for pure Alfvenic waves. We need to analyze the 
relative position of three vectors: magnetic field vector B, wave vector k, and the displacement 
velocity vector v. In what follows, we shall study how the angle 7 between v and B changes with 
plasma (3 . 

It is shown in Alfven &; Falthmmar (1963) that the angle * between v and k can be expressed 
as follows: 

sin cos ,^ , 

where 6 is the angle between k and B, and the phase velocity Vp, is related to the Alfvenic velocity 
Va and the sound velocity Cg through the dispersion relation 

4 - (^i + C's H + CM cos2 9 = 0. (D2) 
Solving this equation for e = vj/vj^, 



m = ^ (1 + /5/2 ± V^(l-/3/2)2 + 2/3sin2ej , (D3) 

where '+' gives the result for fast mode and '— ' represents slow mode. Thus the angle 7 can be 
calculated as 

„ sin 9 cos 9 ,^ ^ . 

7 = 6*- arctan — ^ — (D4) 

cos^ y — e(4) 

and the corresponding plot is shown in Fig. 9. It is evident that for low (3 plasma the velocity v of 
the fast mode is directed nearly perpendicular to B whatever the direction of k, while the velocity 
of the slow mode is nearly parallel to the magnetic field. So the parallel motions we got from slow 
mode are essentially correct, while the perpendicular motions are also subjected to fast mode. A 
more general discussion of the issue is given in CLV02. 
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Fig. 9. — The range of angles between B and v. The dashed area refers to the fast modes, the 
meshed area represents the range where slow modes fall in. The ranges are produced when the 
angle between k and B changes from to tt. 




